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INTRODUCTION 



The J.H. Campbell Plant is comprised of three coal-fired 
operational units. While Units 1 and 2 have used and will continue to 
use Lake Michigan water drawn through Pigeon Lake for cooling purposes, 
Unit 3 draws cooling water from an intake structure located in 
approximately 11 m of water (1.1 km offshore) in Lake Michigan. 
Beginning in September 1980, heated water from all three units was 
discharged through the offshore discharge structure located in 
approximately 6 m of water (0.3 km offshore). In the initial design of 
the benthos survey, 1980 was expected to be an operational year. 
However, due to the late date in 1980 when discharging of warm water 
offshore began, 1980 was considered a preoperational year. As such, 
1980 will be evaluated along with the previous 2 preoperational years, 
1978 and 1979, when comparable data were collected. 

Concurrent with reports on larval, juvenile, and adult fish (Jude 
et al. 1978, 1979, 1980, and 1981), this is the fourth in a series of 
five reports on the benthos that began in 1977. A summary of the intent 
of each report is given in the third of the series (Winnell and Jude 
1980), but for greater detail for 1977 see Jude et al. (1978), for 1978, 
see Winnell and Jude (1979), and for 1979, see Winnell and Jude (1980). 
While the overall purpose of the study design is to determine whether 
regional density and species composition of the benthos collected from 
preoperational years differ significantly from estimates covering the 
operational period, the purpose of this report is to examine all pre- 
operational data for natural, inherent ecological trends, whether they 



be associated with years, months, depths, or regions. Density, 
community structure, and species composition trends will be examined. 
The general distribution of benthos in southeastern and eastern 
Lake Michigan has been summarized in previous reports (Winnell and Jude 
1979, 1980) and need not be restated (in addition to Winnell and Jude 
see Powers and Robertson 1965; Robertson and Alley 1966; Hlltunen 1967; 
Alley 1968; Mozley and Garcia 1972; Mozley and Alley 1973; Mozley 1974, 
1975; Alley and Mozley 1975; Mozley and Winnell 1975; Mozley and 
Howmiller 1977). With each additional year of data, distribution of the 
benthos has become clearer. Previous investigations near the Campbell 
Plant have indicated some regional differences (Truchan 1970, Consumers 
Power Company 1975). During 1978 and 1979, we noted that while depth 
and temporal factors account for much of the variability In benthic 
distributions, regional differences do exist and need to be monitored. 
Analyses of the 1979, plus 1978 and 1979 data combined, further codified 
our perceptions of depth and temporal variability patterns within each 
year and enhanced our conclusions and concerns regarding regional 
differences. Inclusion of 1980 data, which extends the preoperational 
period to 3 yr, will enable us to better vinderstand annual variability, 
particularly as it relates to regional differences. While within-year 
and depth trends over the 3 yr will be compared and assessed regarding 
their permanency, the most important aspect will be confirmation of 
regional differences. A third year of data will clarify regional trends 
noted during the previous 2 yr . 



METHODS 
BENTHOS AND SEDIMENT SURVEY DESIGN 

The survey was composed of 10 stations located along two transects 
perpend Iciilar to the eastern shoreline of Lake Michigan near the 
J.H. Campbell Power Plant, Ottawa County, Michigan (Fig. 1). Along each 
transect, stations were located at 3-, 6-, 9-, 12-, and 15-m depths. 
The first transect represented the treatment area (inner region) located 
0.16 km north of the offshore intake and discharge structures. The 
second transect represented the reference area (outer region) located 
5.0 km north of the offshore intake and discharge structures. The 1980 
survey design was the same as that employed during 1979. 



BENTHOS AND SEDIMENT COLLECTION AND PROCESSING 

Benthic macroinvertebrate and sediment samples were collected on 16 
April, 16 July, and 20 October 1980. Sixty samples were collected for 
benthos and sediments during each sampling month from the University of 
Michigan's R/V Mysis. During 1980, 180 benthic and sediment samples 
were collected. 

Benthos and sediments were collected using a triplex (three- 
chambered) Ponar grab sampler (Mozley and Chapelsky 1973). Each chamber 

2 
of the Ponar samples 0.0165 m . The contents of one side chamber were 

used to estimate numbers of benthic macroinvertebrates occurring in a 

square meter of lake bottom. Contents from the remaining two chambers 

were emptied into a tub and mixed, and approximately 30 g of sediment 

removed for sediment analysis. Six replicates were collected to 

estimate benthic populations and sediments at a station. 
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The portion of the sample used to estimate benthic 
macroihvertebrates was placed in a "funnel-shaped hopper" (see Mozley 
1975 for details) aboard the R/V Mysis. Benthic samples were washed 
through a 0.2-mm mesh net to concentrate animals and remove excess 
sediment and debris. Concentrated samples were stored in labelled 0.5- 
liter Mason jars and preserved with carbonate- buffered, 4% formaldehyde 
solution. Samples were returned to the Great Lakes Research Division 
Benthos Laboratory for sorting and identification. 

Sorting and initial identification of organisms were performed 
using dissecting microscopes (3-30X). Specimens unidentified at the 
genus/species level were mounted on slides with .Amman's lactophenol 
clearing medium and identified using compound microscopes (40-lOOOX) . 

Initial generic identification of chlronomids was determined using 
an unpublished trial key to the chironomids (A.L. Hamilton and 
O.A. Saether, personal communication. Freshwater Institute, Winnipeg, 
Manitoba, Canada and Zoological Museum and Department of Morphology, 
Systematics and Animal Ecology, University of Bergen, Bergen, Norway). 
Where species were determined for chironomid genera, "cf." refers to 
uncertain larval identification at the species level. Most species 
designations concur with reared specimens from the D.C. Cook Nuclear 
Power Plant, southeastern Lake Michigan (see Mozley 1975), which are 
maintained in the Great Lakes Research Division Benthos Laboratory's 
permanent collection. However, since none of the chironomid larvae from 
the Campbell Plant have been reared, identifications at the species 
level have been assigned the uncertainty designator "cf.". The 
designator "gr." refers to a "group" of species inseparable using early 
instar larval morphology and was associated with the genera Chironomus 



and Paracladopelma . Morphology and taxonomy of other chironomid genera 
and species were determined from the following references: Lenz (1954), 
Roback (1957), Curry (1958), Beck and Beck (1969), Saether (1969, 1971, 
1973, 1975, 1976, and 1977), Hirvenoja (1973), Maschwitz (1975), Jackson 
(1977), and Soponis (1977). 

Naidlds were identified using Hiltunen's key to the naidids (see 
Hiltunen and Klemm 1980). Tubificids were identified using an 
unpublished key to aquatic oligochaetes of the Great Lakes 
(J.K. Hiltunen, personal communication; Great Lakes Fishery Laboratory, 
U.S. Fish and Wildlife Service, Ann Arbor, Michigan). Gastropods and 
pelecypods were identified using a key to the mollusks of the Great 
Lakes (Mackie et al. 1980). 

While aboard the R/V Mysis, sediments were stored in "Whirl-pack" 
bags bearing external labels. Standard mechanical sieving of sediment 
samples was performed at the Great Lakes Research Division Sediment 
Laboratory. Folk, Inman, and moment measure statistics were computed 
for each sample collected. Data were expressed in terms of phi units 
following Krumbein (1938). Upchurch (1969), Coakley and Beal (1972), 
and Seibel et al. (1974) indicated that moment measure statistics were 
the "preferred method for deriving sediment textural parameters." Two 
moment measure statistics, mean grain size and standard deviation of the 
mean grain size, were used in this report. Standard deviation was used 
as a measure of sorting, following Seibel et al. (1974). In addition to 
moment measure statistics, percentage of sediments occurring within any 
given sediment grain size based on units of phi was also included. 
Description of sediment grain sizes followed that of Seibel et 
al. (1974), who adapted theirs from the standard Wentworth scale. 



RESULTS AND DISCUSSION 
TOTAL BENTHOS 

The nearshore Lake Michigan benthic community near J. H. Campbell 
Plant was characterized by 113 taxa from 1977 through 1980 (Table 1). 
Of major taxonomic groups present, the most diverse representatives were 
chironomids (44 taxa), pisidia (19 taxa), and naidids (18 taxa). 
Excluding the pilot study year (1977), 101 taxa were observed during 
subsequent years (Table 2). The number of taxa collected during 1980 
(73) was similar to the number observed in 1978 (74), but lower than in 
1979 (80). However, regardless of year, regions have had similar annual 
species richness, with 87 taxa identified from each region. In 1979 
there was a striking difference between regions in the number of taxa 
found at 3 m. Neither this depth nor any of the remaining depths 
exhibited any appreciable differences during 1980, substantiating the 
suspicion that the difference in 1979 was temporary. With respect to 
depth, the number of taxa increased with increasing depth from a low of 
40 taxa at 3 m to a high of 74 taxa at 12 and 15 m. 

Similar to previous years, chironomids, turbellarians , and naidids 
were dominant at 3 and 6 m during 1980, while Pontoporeia hoyi generally 
was the predominant form at 9-15 m (Table 3). Other taxa, in particular 
tubificids, pisidia, and chironomids, constituted a significant portion 
of the benthic community at 9-15 m, in addition to P^. hoyi . Large 
percent occurrence of chironomids at 12 and 15 m in April 1980 was not 

observed to the extent noted during April 1979. 

-2 
The 1980 annual total benthic density estimate (6648 m ) was very 

-2 
similar to that of 1978 (6594 m ) (Table 4). Average benthic densities 

with respect to specific depths or months were quite similar. Mean 
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TABLE 3. Percent occurrence of major taxonomic groups collected in 1980 at 
3 to 15 m among inner (treatment) and outer (reference) regions in the vicinity 
of the J. H. Campbell Plant, eastern Lake Michigan. Percentages expressed in 
terms of total benthos. 















April 
















Inner 










Outer 






Taxon 


3 m 


6 m 


9 m 


12 m 


15 m 


3 m 


6 m 


9 m 


12 m 


15 m 


Chironomidae 


100.0 




47.4 


9.9 


5.6 


28.6 


33.3 


21.5 


20.4 


10.3 


Naididae 




40.0 


2.1 


1.9 


0.7 


28.6 


41.7 


6.5 


2.7 


0.5 


Tubificidae 




20.0 


4.2 


13.5 


6.4 






10.7 


6.5 


16.6 


Enchytraeidae 








0.3 


0.5 






2.0 


8.5 


9.3 


Stylodrilus heringianus 










0.6 










5.9 


Pisidium 






1.1 


7.1 


6.7 


28.6 




2.8 


17.7 


15.9 


Sphaerium 








0.1 












0.5 


Gastropoda 








1.0 


0.7 






0.8 


1.0 


2.6 


Pontoporeia hoyi 






40,0 


65.3 


78.2 




8.3 


53.7 


37.8 


35.9 


Turbellaria 








0.3 


0.5 






1.1 


4.8 


2.2 


Others 




40.0 


5.2 


0.6 


0.1 
July 


14.3 


16.7 


0.8 


0.7 


0.1 








Inner 










Outer 






Taxon 


3 m 


6 m 


9 m 


12 m 


15 m 


3 m 


6 m 


9 m 


12 m 


15 m 


Chironomidae 


84.3 


24.2 


14.1 


4.1 


2.8 


88.2 


51.4 


24.2 


5.3 


2.7 


Naididae 


15.2 


32.2 


19.1 


7.2 


2.1 


11.3 


29.7 


31.3 


18.0 


0.9 


Tubificidae 




3.3 


3.2 


4.7 


4.0 




6.4 


12.8 


14.7 


7.1 


Enchytraeidae 




1.0 


0.2 


1.6 


0.5 




2.0 


0.1 


1.6 


0.4 


Stylodrilus heringianus 






0.1 


0.1 


2.5 








0.1 


6.2 


Pisidium 




0.4 


0.2 


1.4 


6.0 




0.4 


2.4 


4.0 


8.9 


Sphaerium 










0.4 








0.3 


0.1 


Gastropoda 






0.1 




0.1 






0.4 


1.2 


0.4 


Pontoporeia hoyi 


0.4 


38.6 


62.3 


80.8 


80.7 


0.1 


9.9 


27.9 


54.2 


73.2 


Turbellaria 






0.4 


0.1 


0.6 


0.3 


0.2 


0.1 


0.4 


0.1 


Others 




0.4 


0.1 




0.2 






0.7 


0.2 


0.1 












October 
















Inner 










Outer 






Taxon 


3 m 


6 m 


9 m 


12 m 


15 m 


3 m 


6 m 


9 m 


12 m 


15 m 


Chironomidae 


92.2 


99.4 


23.1 


13.5 


4.5 


69.1 


95.2 


43.7 


9.1 


2.6 


Naididae 






51.2 


13.3 


1.3 




1.2 


40.6 


8.3 


1.1 


Tubificidae 






12.4 


27.3 


9.2 






10.5 


31.7 


23.4 


Enchytraeidae 






1.5 


2.4 


4.3 






0.9 


10.8 


3.7 


Stylodrilus heringianus 










2.2 






0.1 




15.7 


Pisidium 


1.1 




1.2 


13.3 


12.6 




1.2 


0.6 


15.2 


18.6 


Sphaerium 








0.2 


0.3 








0.2 


0.4 


Gastropoda 






1.0 


7.6 


5.3 




1.2 


0.4 


5.0 


3.2 


Pontoporeia hoyi 


4.4 


0.6 


6.7 


19.4 


59.0 






1.1 


16.5 


29.9 


Turbellaria 






2.7 


2.8 


0.9 


30.1 


1.2 


1.3 


3.1 


1.2 


Others 


2.2 








0.3 


0.8 




0.9 


0.2 


0.1 
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benthic abundance generally increased with depth during all 3 yr 
(Fig. 2); lowest monthly density occurred in April, and the highest was 
in July, while densities decreased somewhat by October (Fig. 3). While 
several of the regional mean total benthic density comparisons were 
considerably different in 1979, similar estimates in 1980 suggested a 
pattern like that observed during 1978, with no striking differences 
between inner and outer total benthic abundance (Fig. 4 and Appendix I). 
Nonetheless, based on these 3 yr, there was a consistent regional trend 
shown by certain groups comprising the total benthic estimate. The 
group which exhibited the most marked regional trend was P^. hoyi . 

Averaged over all samples in each region (n = 270), P^. hoyi inner region 

-2 
mean abundance was 2808 m , comprising 42% of the benthos; this 

percentage has steadily increased from 1978 (30%) to 1980 (58%). In the 

-2 
outer region, mean abundance of P^. hoyi (2152 m ) was not considerably 

different from the inner region, but percent occurrence was only 29%, 

demonstrating an annual change from 20% in 1978 to 32% in 1980 (Table 

5). Although none of the remaining major taxonomic groups in the outer 

region accounted for large proportions of the difference between 

regions, all were equally abundant or slightly more nijmerous than their 

counterparts in the inner region. Of these groups oligochaetes, 

turbellarians , and chironomids tended to be more abundant in the outer 

region than in the inner region. Most notable was the steady decline in 

percent occurrence of tublficids in the inner region, while in the outer 

region this percentage remained fairly stable from 1978 to 1980. These 

data represent the preoperational status of the regional benthic 

communities, the significance of which will be addressed later. 
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_2 
TABLE 5. Annual percentage of total benthos and mean density (no. m ) for 

major taxonomic groups collected from 1978 to^ 1980^ from the inner (treatment) 

anH n^^1■or■ (-rofa-ra-nno^ Tocr-ionc fn s QD fao-fn-n -.i^ ^ -nam- -f-lio T 11 r"^™-r^V«1 1 



and outer (reference) regions (n 
Plant, eastern Lake Michigan 



90 region yr ) near the J. H. Campbell 













Inner 


Region 










Percentage of 


total 


benthos 
1978- 




Mean density 












1978- 


Taxon 


1978 


1979 


1980 


1980 


1978 


1979 


1980 


1980 


Pontoporeia hoyi 


29.5 


38.8 


58.2 


42.2 


1851.7 


2929.7 


3642.1 


2807.8 


Mysis relicta 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Chiron omldae 


34.6 


12.9 


16.7 


21.4 


2178.2 


974.3 


1047.7 


1400.1 


Naldidae 


14.8 


28.9 


9.5 


17.7 


929.2 


2182.9 


592.5 


1234.9 


Tubificidae 


12.2 


7.3 


6.6 


8.7 


771.6 


549.4 


414.8 


578.6 


Enchytraeldae 


0.9 


1.7 


1.1 


1.2 


59.3 


128.6 


71.4 


86.4 


Stylodrilus herlngianus 


1.4 


0.5 


0.9 


0.9 


85,5 


34.3 


53.9 


57.9 


Oligochaeta 


29.3 


38.4 


18.1 


28.6 


1845.6 


2895.3 


1132.5 


1957.8 


Gastropoda 


1.0 


1.1 


1.2 


la 


62.6 


80.8 


72.0 


71.8 


Pisidium 


3.7 


5.3 


4.8 


4.6 


231.6 


402.7 


300.0 


311.4 


Other Pelecypoda 


0.1 


<0.1 


0.1 


0.1 


3.4 


3.4 


8.8 


5.2 


Hirudinea 


<0.1 


0.0 


0.1 


<0.1 


2.0 


0.0 


4.7 


2.2 


Hydracarina 


0.1 


0.1 


<0.1 


0.1 


8.8 


5.4 


2.7 


5.6 


Hydra 


<0.i 


0.2 


<0.1 


0.1 


2.7 


14.1 


2.0 


6.3 


Turbellaria 


1.6 


3.3 


0.6 


1.8 


100.3 


250.5 


37.0 


129.3 


Other Insecta 


0.0 


0.0 


0.1 


<0.1 


0.0 


0.0 


6.7 


2.2 


Total benthos 










6286.9 


7556.2 


6256.2 


6699.8 










Outer 


Region 










Percentage of 


total 


benthos 
1978- 




Mean density 












1978- 


Taxon 


1978 


1979 


1980 


1980 


1978 


1979 


1980 


1980 


Pontoporeia hoyi 


19.6 


34.1 


32.2 


29.0 


1352.1 


2836.8 


2267.1 


2152.0 


Mysis relicta 


0.0 


<0.1 


0.0 


<0.1 


0.0 


1.3 


0.0 


0.4 


Chironomidae 


32.9 


17.2 


27.0 


25.2 


2269.8 


1432.9 


1898.8 


1867.2 


Naididae 


19.2 


16.6 


11.7 


15.9 


1325.8 


1382.4 


822.8 


1177.0 


Tubificidae 


14.5 


11.1 


U.9 


12.4 


1000.6 


920.5 


839.0 


920.0 


Enchytraeidae 


1.1 


3.0 


2.6 


2.3 


77.4 


249.8 


180.5 


169.2 


Stylodrilus herlngianus . 


2.6 


1.9 


3.7 


2.7 


176.4 


156.9 


263.3 


198.9 


Oligochaeta 


37.4 


32.6 


29.9 


33.2 


2580.2 


2709.6 


2105.6 


2465.1 


Gastropoda 


1.2 


1.3 


1.3 


1.2 


79.5 


106.4 


89.6 


91.8 


Pisidium 


4.4 


6.6 


7.6 


6.2 


301.0 


553.5 


531.9 


462.1 


Other Pelecypoda 


0.1 


0.1 


0.2 


0.1 


9.4 


4.7 


11.4 


8.5 


Hirudinea 


0.1 


<0.1 


0.1 


0.1 


6.7 


3.4 


4.7 


4.9 


Hydracarina 


0.1 


<0.1 


<0.1 


0.1 


8.8 


3.4 


3.4 


5.2 


Hydra 


0.1 


0.0 


0.1 


0.1 


6.7 


0.0 


6.7 


4.5 


Turbellaria 


4.1 


8.1 


1.6 


4.8 


286.2 


671.3 


115.8 


357.8 


Other Insecta 


<0.1 


<0.1 


0.1 


<0.1 


0.7 


0.7 


5.4 


2.3 


Total benthos 










6901.1 


8324.0 


7040.4 


7421.8 
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CHIRONOMIDAE 

While 44 chironomld taxa were identified from the survey area from 
1977 to 1980, only 40 of these taxa were observed during the latter 3 yr 
(Table 1). From 1978 to 1980 the inner region was represented by 35 
chironomid taxa and the outer region by 32 chironomid taxa. Both 
regions were represented by 26 chironomid taxa during 1980 (Table 2). 
The only significant change in chironomid taxonomy from previous years 
has been to tentatively speciate the Polypedilum fallax-gr. individxials 
to Polypedilum cf. tuberculum (see Maschwitz 1975). 

Although dominance of Saetheria cf . tylus and Paracladopelma 
camptolabis- gr . in the survey may occasionally vary, these two 
chironomid taxa were consistently the dominant chironomids (24% and 16%, 
respectively) from 1978 to 1980. For example, in 1980 large densities 
of Psectrocladius cf . s imulans were encountered, making it the most 
numerous chironomid in 1980, although it did not occur very frequently. 
In the inner region, S^. cf . tylus was the dominant chironomid during all 
3 yr, while in the outer region dominance was shared by three different 
chironomid taxa, with S^. cf . tylus the second- or third-most abundant 
chironomid taxon (Appendix 2). 

Although chironomids were represented by a similar frequency in 
samples collected during 1980 (88%) when compared with 1978 (95%) and 
1979 (94%) (Table 4), abundance continued to exhibit considerable 

monthly and depth variability on an annual basis. Mean chironomid 

-2 
density in 1980 (1473 m ) was similar to that observed during 1979 

-2 -2 

(1204 m ); both were considerably less than 1978 levels (2224 m ), 

particularly from 6 to 12 m. Notably different in 1980 was the low 

chironomid density in April and high chironomid density in July 

21 



(Fig. 5), and at the 3-m depth, regardless of month (Fig. 6). Low 
abundance in April was due primarily to changes in S^. cf . tylus density 
resulting from weathering and survival from the previous year. 
Unusually high chironomid density during July and at the 3-m depth were 

produced by two samples which retained second instar Psectrocladius 

-2 
cf . simulans (15,271 and 34,542 m ) in very large abundances and 

accounted for approximately one-third of these density estimates. 

Abundances of this magnitude for P^. cf . simulans are believed to be the 

greatest noted in the Great Lakes and were likely a chance occurrence 

associated with undispersed larvae from a recently deposited egg 

mass(es) . 

Regional comparisons of chironomid abundance in 1980 suggested 

similar regional fluctuations with the exception of October at 9 m 

(Fig. 7). The inner/ outer regional difference at 9 m in October was due 

to greater abundance of Paracladopelma camptolabis- gr . in the outer 

region. Regional densities during 1980 more closely resembled those of 

1979 than 1978. Large differences In chironomid abundance at 15 m in 

April 1980 when compared with 1978 and 1979 were mostly due to decreased 

density of Saetheria cf . tylus . We spectilated that comparatively large 

densities of S^. cf . tylus in April 1979 at 15 m were caused by scouring 

of shallower depths due to spring ice break-up and storm activity. 

However, the weathering process as an explanation for observed 

differences may be incidental to the fact that reproduction and survival 

from 1979 to 1980 may have been considerably less than that which 

occurred from 1978 to 1979, resulting in lower densities in April 1980. 
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NAIDIDAE 

From 1977 to 1980, naidlds were represented by 17 taxa; 16 were 
collected during years after the pilot study year (1977) (Table 1). 
Each region has had 14 naidid taxa identified from 1978 to 1980 of which 
10 taxa were present during 1980 (Table 2). With the exception of the 
very large abundances of Vejdovskyella intermedia in 1979, particxilarly 
in the inner region during July, Piguetiella michiganensis was the most 
abundant and frequently occurring naidid in the survey area^ During 
1980, the naidid Nais variabilis was the second-most numerous naidid 
(Appendix 3), even though previously it was rather scarce in the survey 

area. 

-2 
Lowest annual naidid density was observed during 1980 (708 m ) 

when compared with previous years (Table 4). Likewise, percent 

occurrence in samples collected decreased to 65% in 1980 from a high of 

76% in 1978. The 1980 pattern of monthly and depth density estimates 

most clearly resembled that observed during 1978, with the general trend 

being one of declining naidid density from 1978 to 1980 (Figs. 8 and 9). 

Regional naidid density comparisons indicated the naidid population 

was similarly numerous within regions at each depth considered. In 

addition, annual fluctuations of naidids at each depth were similar in 

both regions (Fig. 10). Although large density decreases for naidids 

were observed in July and October 1980 at 12 m, low abundances were 

observed in both regions, thereby suggesting the causes for the decline 

were operating equally across regions. Finally, large regional density 

discrepancies noted at 6 and 9 m in July 1979 and ascribed to 

construction activity impact appeared to have been ameliorated in 1980 

as had been expected . 
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TUBIFICIDAE 

Tubificids were represented by 13 taxa from 1977 to 1980, but only 
11 of these were collected during the latter 3 yr (Tables 1 and 2). 
Seven tubificld taxa were observed in each region during 1980, bringing 
the total number of tubiflcid taxa in the inner region to 9 and the 
outer region to 10. Potamothrix moldaviensis was the most abundant 
tubificid during each year and was similarly numerous during all years 

in each region, except 1979 in the inner region. 

-2 
Annual mean tubificid density (627 m ) and percent occurrence in 

benthlc samples (58%) during 1980 were the lowest annual values observed 

in the survey area from 1978 to 1980 (Table 4). Although there was a 

slight decrease in annual tubificid abundance, the distributional 

pattern of mean monthly densities and mean depth densities continued to 

be very similar during each year (Figs. 11 and 12). 

While annual tubificid abtmdances decreased in both regions from 

1978 to 1980, the decline was comparatively greater in the inner (46%) 

than in the outer (16%) region. Although regional density estimates for 

tubificids indicated greater abundance in the outer as opposed to the 

inner region at 12 and 15 m during most months sampled (Fig. 13, 

Appendix 3), most of these differences were not extreme. However, the 

most disparate of these regional differences occurred in October at 15 m 

where tubificid numbers continued to increase in the outer region and 

decrease in the inner region from 1978 to 1980. Continuation of the 

October, 15-m regional difference can at present be ascribed to no other 

cause than inherent differences between the two regions. Remaining 

regional differences were likely of little significance, given that 

tubificid densities are normally very patchily distributed. 
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ENCHYTRAEIDAE 

Enchytraeids were represented by at least two species in 34% of the 

-2 
samples collected during 1980, having an annual mean density of 125 m 

(Table 4). Both frequency of occurrence and abundance population 

parameters for 1980 enchytraeids were intermediate to those of previous 

years, with 1979 being by far the most abundant year for enchytraeids. 

While deviating slightly from previous years, enchytraeid abundance with 

respect to month and depth generally concurred with distribution 

patterns observed in 1978 and 1979 (Figs. 14 and 15). Regionally, 

greatest enchytraeid density differences during 1980 were at 15 m in 

April and 12 m in October. Both of the regional density differences 

were of temporary duration. The former difference was no longer evident 

by July 1980. The latter regional difference, although requiring April 

1981 data for confirmation, was suspected to follow a pattern similar to 

that noted from October 1979 to April 1980 at 12 m, i.e., minimal 

difference, followed by no difference in July 1980 at 12 m (Fig. 16, 

Appendix 3) . 



STYLODRILUS HERINGIANUS 

Occurrence of Stylodrilus heringianus in the survey area was 
limited almost exclusively to 15 m in 1980 as well as during other years 
(Appendix 3). S^. heringianus was collected in 19% of the samples during 
1980, which was comparable to values observed in previous years (Table 

4). Greatest annual mean abundance of ^. heringianus was observed 

-2 
during 1980 (159 m ), although this value was not much greater than 

-2 -2 
those observed during 1978 and 1979 (131 m and 94 m , respectively) . 

The same distributional pattern was observed in the survey area with 
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respect to month (gradual Increase In abundance from April to October In 
response to reproductive activity) and depth (greatest abundance at 15 
m) during all 3 yr (Figs. 17 and 18). However, based on a 3-yr 
estimate, S^. heringianus was consistently more numerous in the outer 
(199 m ) than in the inner (58 m ) region (Table 5). 

Comparison of S^. heringianus regional mean abundance suggested no 
disparate trends during 1980, except at 15 m in October (Fig. 19). S^. 

heringianus density increased over two- fold at 15 m during October in 

-2 -2 

the outer region from approximately 1000 m in 1978 to nearly 2400 m 

in 1980. In the inner region, however, abundance of S^. heringianus 

-2 -2 

decreased from approximately 1000 m in 1978 to less than 250 m 

during 1980. In contrast to October, regional S^. heringianus densities 

during April and July at 15 m were similar. Therefore, even though 

October regional differences at 15 m were increasingly disparate with 

each successive year, these dissimilarities have not resulted in any 

consistent regional trends during other months and years. No 

explanation can be offered at present to clarify the cause of this trend 

other than inherent differences between regions. 



TURBELLARIA 

Three unknown tur.bellarian forms were identified from the survey 
area. As in 1979, the species 1 form was dominant at 3 and 6 m, with 
the species 2 form generally most numerous at remaining depths (Appendix 
1). Species 3 was also found at 9-15 m but in very low densities. 
Aside from these similarities, dissimilarity appeared to be the 
standard, as annual variation of turbellarian density was very extreme 
from 1978 to 1980. Frequency occurrence decreased from a maximinn 
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observed in 1979 (64%) to 34% in 1980, almost the same as that observed 
during 1978. Average annual abundance decreased substantially from 461 
m~ in 1979 and 193 m~^ in 1978 to 76 m~^ in 1980 (Table 4). Although 
October abundance during 1980 was lower, the turbellarian seasonal 
distribution pattern in 1980 was similar to 1978 (Fig. 20). Turbellarian 
densities ^across depths were very similar during 1980, with fewest 
collected at 6 m. While resembling the 1978 density trend with depth 
more closely than that of 1979, the number of turbellarians occurring at 
any given depth was unpredictable (Fig. 21), except at 6 m where lowest 
abundances occurred during each year. 

Outer region turbellarian densities (358 m ) were consistently 

-2 
greater than inner region values (129 m ) based on a 1978-1980 average 

(Table 5). While this regional difference persisted for turbellarians 

in general in the survey area, it was most evident at 3 m during July 

and October during all 3 yr (Fig. 22). Other depth and month 

combinations demonstrated temporary regional disparities which were 

minimized during succeeding months or years. 



PELECYPODA 

Three genera and 17 species of pelecypods were found in the survey 
area from 1977 to 1980. The majority of pelecypod species belong to the 
genus Pisidium (14 species), with two species of Sphaerium and a single 
species of Musculium also represented (Table 1). Although occurring in 
54% of samples collected, pelecypods comprised only a small portion of 
the benthic population (5.6%) (Table 4). Pisidium was the dominant 
pelecypod genus, comprising 98.2% of the pelecypod abundance or 5.5% of 
the benthos (Table 4). 
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Plsidlum 

A shift in species dominance from P^. nitidum nitidtim in 1978 to 
P^. casertanum in 1979 was more pronounced during 1980, although rank 
ordering of the five dominant pisidia was the same as that observed 
during 1979. Percent occurrence of P^. n. nitidum continued to decrease 
from 1978 (31%) to 1980 (13%), and P. casertanum increased from 1978 
(23%) to 1980 (42%). Percent occurrence of P. fallax, which increased 
from 1978 to 1979, did not change much from 1979 (19%) to 1980 (15%). 
The change in species dominance occurred similarly in both regions 
(Appendix 4) . 

Distributional patterns of Pisidium followed the same general 
patterns established during the previous 2 yr. Mean monthly pisidia 
abundance increased from April through October during 1980 (Fig. 23). 
Increased pisidia density from April to October was mostly due to annual 
reproduction, as October specimens were primarily small individuals. 
Average density of Pisidium at each depth in 1980 was most similar to 
that observed during 1978, but nonetheless increased with increasing 
depth during all 3 yr. Very few pisidia were observed at 3 and 6 m in 
any year (Fig. 24), likely due to coarser substrates and higher degree 
of turbulence and shifting sands not conducive to the filter feeding 
activity of clams. 

Regional comparisons of mean Pisidium abundance indicated sporadic 
differences, although the general trend indicated densities were 
consistently less in the inner when compared with the outer region 
(Fig. 25). Since these discrepancies were minimal in either succeeding 
months or years at a given depth, and since regional species composition 
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was very similar, little significance can be assigned to observed 
regional differences at this time. 

Sphaerium and Musculium 

The sphaeriid species Musculium transversum , previously assigned to 
the genus Sphaerium , has been reassigned to Musculium (Mackie et 
al. 1980). Consequently, there are now two genera where earlier reports 
listed only Sphaerium . M. transversum was observed only during 1979. 

S^. striatinum has continued to be the most abundant Sphaerium 

-2 
species in the survey area (6m ). When all species of Sphaerium and 

Musculixim were combined to represent sphaeriids exclusive of Pisidium , 

similar annual densities were observed from 1978 to 1980, with maximum 

yearly abundances occurring at 12 and 15 m (Fig. 26). 

Sphaerium nitidum continued to be poorly represented at depths 

sampled during the survey years 1978 to 1980. Annual densities of 

-2 
Sphaerium nitidum have never exceeded 1.3 m . Since 1978, S^. nitidum 

density has decreased until in 1980 no specimens were collected 

(Appendix 4). However, it is likely that 12 to 15 m is the nearshore 

extent of a S^. nitidum population which is more successful in deeper 

waters. During 1977 at 20 and 25 m, S^. nitidum was observed in average 

-2 
densities of approximately 100 m 



GASTROPODA 

As during previous years, Valvata sincera was the most abundant of 
five gastropod taxa in the survey area during 1980 (Appendix 4). 
Lymnaea sp. continued to be the second-most nimxerous gastropod as was 
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the case during 1979. Gastropods were again represented in one-third of 

the samples collected and had an annual mean density during 1980 (81 

-2 
m ) similar to those found in previous years (Table 4). 

Although average gastropod density at each depth in 1980 was 

similar to other years and followed the same distributional pattern 

-2 
(i.e., low numbers at depths less than 12 m and 130-250 m at 12 and 15 

m) , mean monthly density varied (Figs. 27 and 28). Regardless of annual 

or monthly density fluctuations for gastropods, regional density trends 

were similar for nearly all depths within each month sampled (Fig. 29). 



PONTOPOREIA HOYI 

-2 
Annual mean abundance of chironomlds (2224 m ) and oligochaetes 

-2 -2 

(2213 m ) was greater than that of P. hoyi (1602 m ) during 1978. 

However, during subsequent years, P^. hoyl was the dominant organism 

collected in both regions of the survey area, having an anniial mean 

density of 2883 m"^ in 1979 and 2955 m"^ in 1980 (Tables 4 and 5). P. 

hoyl abundance in 1980 most closely approximated that of 1979, with 

monthly and depth distributional patterns resembling those of previous 

years, i.e., high densities in July following release of young in April 

and May, and with respect to depth, relatively few individuals at 3 and 

6 m (Figs. 30 and 31). 

An increasing annual density trend was observed from 1978 to 1980 

In the inner region that was not observed concurrently in the outer 

region. Yearly average abundance of P^. hoyi increased nearly two-fold 

in the inner region from 1978 (1852 m~^) to 1980 (3642 m~^). In 

addition, £. hoyi comprised 58% of the inner region benthlc abundance in 

1980 which compared with only 30% in 1978. 
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In the outer region, following a large increase in P^. hoyi mean 

-2 -2 

density from 1978 (1352 m ) to 1979 (2837 m ), there was a subsequent 

-2 
decrease in abundance during 1980 to 2267 m . However, percent 

contribution of jP. hoyi to the outer region benthic abundance remained 

at approximately 33% during the last 2 yr (Table 5). Consequently, 

while P^. hoyi comprised an average 29% of the benthic community 

abundance in the outer region from 1978 through 1980, the inner region 

benthic community was comprised by 42% P^. hoyi . The significance of 

this regional difference was largely that outer region P^. hoyi 

percentage of the benthos was similar across years while that in the 

inner region increased with each successive year. Whether the 

percentage of jP. hoyi in the inner region benthic community continues to 

rise is yet to be determined, but the inner region has consistently 

supported greater annual mean densities of P^. hoyi when compared to the 

outer region. When averaged over all years the inner region P^. hoyi 

-2 -2 

average density was 2808 m which compares with 2152 m in the outer 

region. 

Greatest differences between regions occurred during July at 9 and 

12 m in 1980 (Fig. 32, Appendix 1). While a regional difference at 9 m 

in July was first noted in 1980, the difference at 12 m in July 

continued a trend first observed during 1979. In addition, at 15 m in 

July there was a reversal of the regional P^. hoyi density trend during 

1979, with the inner region P^. hoyi density becoming greater than that 

of the outer region in 1980. However, in all examples regional 

differences in jP. hoyi density were nearly non-existent in October 

samples. 
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In 1978 and 1979, analysis of P^. hoyi size classes indicated an 
apparent trend which was characterized by individuals of the inner 
region at 9 and 12 m in April and July being somewhat advanced, i.e., 
one size class larger (see key in Fig. 33 for size-class definition), 
than those of the outer region. In both years this difference was 
minimized by October (see Winnell and Jude 1979, 1980). A similar 
analysis in 1980 indicated very little regional difference in size class 
attained by F_. hoyi at given depths for specific months (Fig. 33). 
Regional size-class differences were most evident in 1978, and although 
demonstrable in 1979, relative to 1978 they were less pronounced, and 
appeared to have become negligible during 1980. 

There continued to be an onshore to offshore size-class trend 
associated with depth in the range 9 to 15 m. Most P^. hoyi individuals 
collected at 9 and 12 m in April and July were one size class larger 
than those at 15 m. However, these differences were non-existent by 
October. This depth trend occurred in both regions. 



SUBSTRATE DISTRIBUTION 

Based on yearly mean percentage of sediment parameters, there has 
been virtually no difference in composition of the substrate from 1978 
to 1980 (Table 6). The average substrate encountered in each year was 
well to moderately sorted fine sand, with 93% of substrate types lying 
in the medium to very fine sand range (1-4 units of phi). 

Within each year, when comparing regional mean percentages of each 
sediment parameter at specific depths within months, some differences 
were noted. The greatest difference observed occurred during July 1979 
in the inner region where a large increase in silt and clay substrates 
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Fig. 33. Percent distribution of P^. hoyi size classes in the inner 
and outer regions at 9-15 m during April, July, and October 1980. 
Samples were collected from eastern Lake Michigan near the J. H. 
Campbell Plant. (* = <2%) . 
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was observed compared with previous sampling periods. As this regional 
substrate difference was negligible by October 1979, little regional 
difference was expected during 1980. Comparison of regional substrate 
trends during 1980 confirmed this view as few differences were present 
(Table 7). Most regional differences were minor percent fluctuations 
between medium- and fine-sand categories, although during July the inner 
region did have nearly twice the percentage (averaging 13.2%) of very 
fine sand (3-4 phi), than the outer region (6.7%). Again in October at 
3 m there was a considerable regional difference, iwith the inner region 
comprised of a finer substrate (mean phi » 2.89) than the outer region 
(mean phi » 1.89). However, since medium to very fine sands were the 
predominant substrates expected to be encountered in the survey area in 
proportions that may be temporally and spatially variable, there was 
little reason to regard regional differences obseri^ed during 1980 with 
the significance of that ascribed to July 1979. Consequently, regional 
substrate types during 1980 were not thought to represent any definitive 
deviant trends from the expected sediment distribution patterns 
characteristic of previous years. 

SUMMARY AND PROJECTED ANALYTICAL AND INTERPRETI^TE PROBLEMS IN 1981 



One purpose of this 4-yr study is to assess the effect of heated 
waste water (discharged from cooling processes in the Campbell Plant) on 
the macrobenthic community. At this juncture in data collection and 
analysis, two problems will complicate evaluation of thermal impact. 
The first problem is directly related to the conceptual design of the 
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TABLE 7. Average percent composition of sediments distributed among 
sediment grain sizes, average mean grain size (x) , and average standard 
deviation of mean grain size (SD) for sediments collected in 1980 at 
3 to 15 m in the inner (treatment) and outer (reference) regions (n = 6) 
near the J. H. Campbell Plant, eastern Lake Michigan. 
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survey's analytical methodology. Wlnnell and Jude (1980) defined the 
application of Johnston's (1973, 1974) mixed-model, nested analysis of 
variance (ANOVA) design to the present data for estimation of heat 
effect and precision of detection. The ANOVA was established as a 
balanced design with five essential factors; construction time (before 
and after offshore theinnal discharge), year (nested within construction 
time), month, depth, and region [treatment (inner) and reference 
(outer)]. However, the offshore thermal discharge of Unit 3 did not 
commence until September 1980, instead of prior to April 1980 as had 
been anticipated. Consequently, the balanced design was altered to an 
unbalanced design with 3 yr of preoperational data (1978 through 1980) 
and 1 yr of operational data (1981) from the previous assumption of 2 
preoperational (1978 and 1979) and 2 operational (1980 and 1981) years. 
Therefore, it will be necessary to evaluate the practicability and 
propriety of the resultant unbalanced design's application to assessment 
of heat effect. Given the conceptual change in the design, other models 
wj.ll be considered that may be more applicable. 

A second problem that may limit or at least influence 
interpretation of significant interactions of construction time and 
region (heat effect) is the presence of preoperational regional 
differences for certain components of the benthic community. Analysis 
of 3 yr of annual mean density data gathered from 1978 to 1980 strongly 
suggests that P^. hoyi is becoming steadily more numerous and is 
comprising increasingly greater proportions of the benthos in the inner 
region when compared with the outer region. Several other major 
taxonomic groups, e.g., most oligochaetes, turbellarians, and 
chironomids, are somewhat more numerous in the outer region when 

87 



compared with the inner region, although the size of these differences 
was not as large as those noted for P^. hoyi « Several of the benthic 
community components, i.e., tubificids, S^. heringianus , and Pisidivim , 
had consistent regional density differences during the last 2 yr at 15 m 
in October. Although recruitment from reproductive activity produce a 
large proportion of individuals encountered during October, whether 
regional differences were coincidental or causative remains 
undetermined. As these general trends appeared to be operating in the 
survey area, specific levels of comparison, (i.e., regional comparisons 
at each depth within each month sampled) are quite variable, with large 
density differences characteristic of 1 mo minimized in ensuing months 
and/or years. Annual abundance differences between regions for some 
taxa may or may not be adequately expressed in specific regional 
comparisons. Therefore, until the ANOVA is completed, no conclusive 
statements can be made regarding the significance of trends. However, 
presence of regional trends in preoperational data will require careful 
evaluation of ANOVA results. 

Sedimentary trends observed in the 3-yr period from 1978 through 
1980 generally lack annual distinctions. Substrates at 3-15 m were well 
to moderately sorted fine sand, with over 90% of the substrates being 
medium to very fine sand. While there was an expected depth-dependent 
distribution of sediments, overall trends have failed to indicate any 
regional differences in substrates, excepting the unusually high amounts 
of silt and clay observed in the inner region during July 1979 
construction activity. 
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